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Abstract 

By considering a first-order variation in electroaccepting and electrodonating powers, or, induced by a change 
from gas to aqueous solution phase, the solvent effect on or for charged ions is examined. The expression of 
electroaccepting and electrodonating powers in the solution phase, ar s , is obtained through establishing the 
quantitative relationship between the change of the or due to the solvation and the hydration free energy. It is 
shown that cations are poorer electron acceptors and anions are poorer electron donors in solution compared to 
those in gas phase. We have proven that the scaled aqueous electroaccepting power, co + s , of cations can act as a 
good descriptor of the reduction reaction, which is expected to be applied in the design of solution reactions. 
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Introduction 

With the rapid development of functional materials, 
novel micro/nanostructures of the materials are highly 
demanded to obtain advanced properties, which can be 
achieved by the rational design of solution-phase chemi- 
cal reactions [1-6]. Therefore, it is of significance to 
thoroughly understand the reactivity of chemical species 
and the mechanism of chemical reactions to further rea- 
lize the solution reaction design. Among many chemical 
reactivity indices, one quantity of importance is the elec- 
trophilicity, oj, introduced by Parr and co-workers [7]. 
They defined co as 



(1) 



where \i is the chemical potential and T) is the chemical 
hardness of an TV-electron system with total energy, E, 
defined as fi = (3 2 £/3 2 TV) v(r) and tj = (3 2 £/3 2 TV) v(r) . This 
index has been found to be helpful in analyzing the reac- 
tivity behaviors of a variety of compounds as well as the 
reaction mechanisms of diverse chemical processes [8,9]. 
As an important contribution to the co, Gazauez et al. 
[10] argued that from a chemical perspective, it would 
make sense to differentiate the response of a system to 
the electron acceptance from the electron donation 
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grounded on that the left and right derivatives of the 
total energy, £ DFX (TV), of an TV-electron system with 
respect to the integer electron number, TV, are different. 
By introducing an electron bath of nonzero chemical 
potential, ^bath> with which the chemical species can 
exchange electrons, they proposed electroaccepting [co + ] 
and electrodonating [co ] powers as the following equa- 
tion: 



(M ± ) 2 ~ (Mbath) 2 

2rj± 



(2) 



where the chemical potential, fc + , and the chemical 
hardness, rf, were defined as 



(J + 3A) + I — A 

V- = : — f = r l = —r- 



(31 + A) 
[i = = i] 



I — A 



(3a) 



(3b) 



where I and A are the ionization potential and the 
electron affinity, respectively. A larger value of or corre- 
sponds to a larger capability of accepting charges, 
whereas a smaller value of co' implies a larger capability 
of donating charges. 

Although some chemical phenomena have been ratio- 
nalized by establishing the quantitative structure-reactiv- 
ity relationships using these reactivity indices in the gas 
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phase [11], the presence of solvent is bound to affect the 
reactivity behaviors of chemical substances. Therefore, 
studies on the reactivity indices such as m and co in 
solution are quite necessary to reveal the accurate reac- 
tivity of chemical species in solution and further predict 
and design the solution phase reactions [12]. While sev- 
eral theoretical calculations about the solvent effect on 
the m for various chemical species have been performed 
[12-14], the solvent effect on the w ± which are regarded 
as better descriptors of the donor-acceptor type interac- 
tions [10] has not received much attention to date. In 
this work, the solvent effect on the co* is estimated by 
establishing a linear relationship between the change of 
the m~ due to the solvation and the hydration free 
energy, AGh y d- The values of aqueous electroaccepting 
power, eo + s , of 39 metal cations are quantitatively calcu- 
lated, which are proven to be appropriate descriptors for 
the reduction reactions. 

Method 

In a previous study, Perez et al. [14] examined the sol- 
vent effect on the electrophilicity index, m, by introdu- 
cing a first-order finite variation in the m due to the 
solvation 



A<w e - 



where A^ g _> s and Ar7 g _> s are the variations in ^ and r\ 
from the gas to solution phase, respectively. 
They rearranged the first contribution as 

■ ft) ■ [SS A ^ s 3 AEi ™ = 2AGso,v ® 

where AE ins is the insertion energy of the solute going 
into the solvent which is suggested as twice the solva- 
tion energy. 

The second contribution in Equation 4 is rewritten as 



2V>1 E / ' s iVoJUgAW 1 VlgAN/ \l) g ") 



1 / ft. \ AN„ malI 



Finally, they deduced the expression of Am 



AG SC 



(6) 



(7) 



where Aco g ^ s showed a linear dependence on the sol- 
vation energy, AG so i v , with a regression slope, y. They 
used 18 well-known electrophilic ligands including hard 
electrophiles such as Li + and Na + to test this linear cor- 
relation and obtained good results (R = 0.9925, y = 
1.00765 at B3LYP/6-311G"and R = 0.9918, y = 0.96843 
at HF/6-311G**levels of theory). 



Herein, we reconstruct the second contribution in 
Equation 4 which will directly lead to a quantitative 
expression for Aco g ^> s with a definite slope value, y. 



=4(--) 2 [(^)-(^)]-^-) 2 

1, .2 (0 — Ms) — (0 — Ito) 1, ,1^-^ 
, AN„J — ^ (AN ema x) — — 



(8) 



Substitution of Equations 5 and 8 into Equation 4 
leads to the expression of Aca g _^ s . 

Afflg^ s = Att> ljg ^ s + Ati> 2 ,g^ s = 2AG so iv - AG so i v = AG so i v (9) 

Therefore, the global electrophilicity, m s , in solution 
can be calculated by 



a), = u>a + Aa>„ 



aj g + AGsoiv. 



(10) 



It should be noted that one key assumption in our 
approach is A/V s , max ~ A7V g> max which could be justified 
by the data of Table 1 in Perez's work [14]. Our result, 
7=1, has turned out to be fairly consistent with Perez's 
regression value, i.e., y = 1.00765 and y = 0.96843, 
which thus approve the reasonableness of our approach 
to dealing with Am 2 , g _> s - 

Further, we try to extend our approach to examine the 
solvent effect on the or. For the charged ions, we sup- 
pose that the chemical potential, ^ ± bath> of the electron 
bath equals that of the parent atoms of ions since the 
charged ions become neutral atoms after accepting or 
donating the maximum amount of electrons. In addi- 
tion, as the solvent only has little effect on the chemical 
potential, fi, of the neutral species [12,13,15], there exists 
a relationship as ^ ± bat h = ^s, atom » f* s , atom- The ion 
exchanges electrons from the bath to the point that its 
chemical potential, , equals the value ^ ± b a th with the 
maximum amount of electron flow: 

AN ± max = ( M ± bath -M ± )/^ ± . (11) 

The first-order variation in the of leads to the follow- 
ing equation: 



(12) 



The first part of Equation 12 in terms of the variation 
in pr is given by 





r A£ g i 


AN± S 






L^sJ 



A/iV, = AChyd (13) 



where the energy change Aa) ± l g _ >s due to the variation 
of the chemical potential from the gas to solution phase 
can be represented by AGh y d [16,17]. 

The second part of Equation 12 in terms of the varia- 
tion in rf is given by 
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Table 1 Calculated electroaccepting power, o>* s , 


in aqueous solution and the absolute reduction potential, £° abs 


M z+ 






„+ b 
7] 


AN*,. max 

y, nicJX 


AG hyd 


Ara + 


C0+ 5 




Li + 


-4 044 


-22.954 


35.1 24 


0.538 


-481 


-2.053 


5.214 


1 .223 


Na + 


-3.854 


-15.677 


21.075 


0.561 


-375 


-1 .466 


4.013 


1.550 


K 


-3.255 


-1 1 .163 


1 3.645 


0.580 


-304 


-1.1 16 


3.062 


1.339 


Rb + 


-3.133 


-9.953 


1 1.551 


0.590 


-281 


-0.998 


2.865 


1.339 


Cs + 


-2.920 


-9.1 95 


1 0.603 


0.592 


-258 


-0.912 


2.673 


1.340 


Aa + 


-5.682 


-1 1 .055 


6.957 


0.772 


-440 


-1.108 


5.354 


5.062 


Cu + 


-5.795 


-10.867 


6.282 


0.807 


-535 


-1 .294 


5.433 


4.783 


T! 


-4.581 


-9.688 


7.160 


0.713 


-310 


-0.847 


4.242 


3.923 


In + 


-4.340 


-9.057 


6.542 


0.721 


-296 


-0.799 


4.032 


4.123 


Be 2+ 


-6.992 


-52.133 


67.843 


0.665 


-2404 


-10.787 


8.883 


4.546 


Mg 2+ 


-5.735 


-31.313 


32.554 


0.786 


-1838 


-7.780 


6.774 


3.814 


Ca 2+ 


-4.585 


-21.632 


19.521 


0.873 


-1515 


-6.187 


5.260 


2.846 


St 2 * 


-4.271 


-18.995 


15.930 


0.924 


-1386 


-5.567 


5.185 


2.746 


Ba 2+ 


-3.909 


-16.378 


1 2.748 


0.978 


-1 258 


-4.963 


4.958 


2.686 


V 2+ 


-5.060 


-18.291 


7.347 


1.801 


-1825 


-6.841 


14.187 


6.266 


Cr 2+ 


-5.075 


-20.105 


7.237 


2.077 


-1860 


-7.206 


1 8.940 


6.726 


Mn 2+ 


-5.576 


-20.147 


9.014 


1 .61 7 


-1 770 


-6.634 


14.157 


6.1 86 


Fe 2+ 


-5.927 


-19.804 


7.232 


1 .91 9 


-1 848 


-6.71 1 


1 7.976 


7.646 


Co 2+ 


-5.91 1 


-21.188 


8.208 


1.861 


-1922 


-7.181 


18.037 


7.972 


Ni 2+ 


-5.730 


-22.424 


8.51 1 


1 .962 


-1 998 


-7.708 


1 9.906 


8.012 


Cu 2+ 


-5.795 


-24.429 


8.275 


2.252 


-2016 


-7.969 


26.064 


9.206 


Zn 2+ 


-7.046 


-23.405 


10.879 


1.504 


-1963 


-7.1 10 


1 5.784 


7.001 


Cd 2+ 


-6.745 


-22.051 


10.286 


1.488 


-1 736 


-6.244 


15.180 


7.720 


Ha 2+ 


-7.828 


-22.618 


7.722 


1.915 


-1766 


-5.984 


23.174 


10.234 


Sn 2+ 


-5.508 


-18.600 


7.936 


1.650 


-1496 


-5.457 


1 4.429 


8.251 


Pb 2+ 


-5.562 


-19.258 


8.453 


1.620 


-1434 


-5.285 


1 4.824 


8.274 


Pd 2+ 


-6.253 


-22.805 


6.750 


2.452 


-1920 


-7.222 


28.406 


10.356 


Sm 2+ 


-4.233 


-14.153 


6.165 


1.609 


-1375 


-4.994 


9.797 


3.186 


Eu 2+ 


-4.253 


-14.668 


6.835 


1.524 


-1391 


-5.1 18 


9.296 


2.926 


Yb 2+ 


-4.691 


-15.395 


6.437 


1.663 


-1510 


-5.441 


1 1 .259 


2.926 


Al 3+ 


-4.489 


-51.334 


45.772 


1.023 


-4531 


-21.427 


7.139 


7.761 


Ga 3+ 


-4.499 


-39.033 


16.645 


2.075 


-4521 


-20.728 


24.430 


1 1 .202 


ln 3+ 


-4.340 


-34.523 


1 2.985 


2.324 


-3989 


-18.073 


27.093 


11.775 


Sc 3+ 


-4.921 


-36.940 


24.366 


1.314 


-3801 


-17.073 


10.431 


6.699 


Y 3+ 


-4.663 


-30.539 


20.039 


1.291 


-3457 


-15.179 


7.550 


5.679 


La 3+ 


-4.183 


-26.870 


15.386 


1.475 


-3155 


-13.805 


9.090 


5.649 


Fe 3+ 


-5.927 


-36.689 


1 2.074 


2.548 


-4271 


-18.557 


35.731 


1 2.669 


Co 3+ 


-5.91 1 


-37.950 


8.900 


3.600 


-4503 


-19.701 


59.247 


14.148 


Au 3+ 


-6.919 


-41.575 


8.350 


4.150 


-4416 


-19.075 


81.561 


1 7.349 



3 A t+ bath = M atom = — (3/m + j 4 ^ — 3/m / 4 because of the neglectable values of electron affinities, A M , for atoms. b jj + and rj + are obtained 
from Equation 3a. c The product of the charge number, z, and the absolute reduction potential, E° abs . The values of P abs are calculated form Equation 18 where 
the absolute standard hydrogen electrode potential, P S he = ' s 4.263 eV according to Marcus [21]. 
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M + bath = I 1 atom = -(3/ M + A M )/4 = -3J M /4 (14) 



Combining Equation 13 with Equation 14 yields 



A&) ± g_* s = Aw ± i j g_ >s 4- Aoj ± 2,g 



/ _ 1 /^bath + ^ 
\ 2 f±. 



, 1 M bath + M g \ A ± 
1 - i I A(B i, g . 



2 M ± 



(15) 



AG hy(i . 



Therefore, the electroaccepting and electrodonating 
powers in solution, <±r s , can be calculated by 



; ^rV - (16) 
1 I- 1 8 / 



Results and discussion 

According to Pearson's viewpoint that cations are elec- 
tron acceptors and anions are electron donors [16], we 
pay attention to the oj + s for cations and m' s for anions. 
By using Equation 16, the m + s values for 39 metal 
cations with charges from +le to +3e are calculated and 
summarized in Table 1. From Table 1, we find that the 
solvation weakens the capacity of cations to accept elec- 
trons due to the negative values of Aco + , in agreement 
with the previous conclusions [12-16]. Unfortunately, it 
is impossible to quantitatively calculate the of s values 
for anions so far due to the absence of experimental 
electron affinities needed in Equation 3b. Herein, these 
values can be qualitatively estimated: 



> 0) e + I 1 



1 ft g + IX 

2 n~„ 



_ 1 M bath + M 

2 fi- g 

AGhyd > <«~g- 



AGhyd 



(17) 



Since a larger value of co' s implies a smaller capability 
of donating charges, we can conclude from Equation 17 
that the solvation also weakens the capacity of anions to 
donate electrons, which agrees with the general view- 
points [12-16]. 

Many liquid-phase chemical reactions involve the elec- 
tron-transfer steps, and a key thermodynamic variable that 
describes the tendency of chemical species in solution to 
gain or lose electrons is the redox potential. The quan- 
tum-chemical computation approach to electrochemistry 
has become available very recently [18]. However, the esti- 
mation of redox potential by the quantum-chemical calcu- 
lations is a great challenge due to the complexity of the 
processes involved in a typical electrochemical reaction 
[19]. For example, the complicated diffusion and adsorp- 
tion processes on the electrode surface which should be 
necessarily taken into account in the quantum-chemical 
modeling of the reduction-oxidation reaction lead to the 
considerable system size and thus require strong comput- 
ing power. Therefore, previous studies mainly focus on the 
one-electron reduction reactions between different oxida- 
tion states of transition metals to avoid modeling of an 
electrode-solution boundary [18-20]. In this work, we try 
to use the m + s to describe the many-electron reduction 
reaction including both main- and sub-group metal 
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Figure 1 Plot of the m* s versus zE° abs . 
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cations. According to the reaction formula M + (aq) + zl 
2H 2 (g) = M (aq) + zH + (aq), the absolute reduction poten- 
tial, £° a bs» can be calculated by 

^bs = £sHE.O+£°abs(SHE) (18) 

where £° abs (SHE) is the conventional reduction poten- 
tial and £°she = o is the absolute standard hydrogen 
electrode potential. Note that the m + s is the energy low- 
ing associated with a maximum amount of electron flow 
between two species; it is reasonable to establish a cor- 
relation between z£° abs and co + s . A good relationship 
shown in Figure 1 approves that our co + s can act as an 
appropriate descriptor of the many-electron energy 
change. Moreover, this method is more simple and con- 
venient compared to the quantum-chemical approach to 
the estimation of the £° a bs- 

Except for the reduction reaction, the £0*5 can also be 
expected to qualitatively and quantitatively predict other 
properties of ions in connection with ligand binding, 
hydrolysis processes, and stability of coordination com- 
pounds, etc. In addition, compilation of experimental 
data on solvation energies in nonaqueous solutions will 
make it possible to evaluate the corresponding electroac- 
cepting and electrodonating powers, w , which will 
undoubtedly lead to the deeper understanding of the che- 
mical reactivity of ions in these media. 

Conclusions 

By reconstructing a first-order variation of the m due to 
the solvation, the linear relationship between the change 
in the m and the solvation energy is reproduced, which 
suggests that our method is theoretically reasonable. The 
solvent effect on the electroaccepting and electrodonat- 
ing powers, or, for charged ions is examined, and a defi- 
nite quantitative expression for the aqueous or s is 
established. It is found that the solvation weakens the 
capability of both electron-accepting power of cations 
and electron- donating power of anions. A good relation- 
ship between the co + s and £° a b s shows the validity of the 
electroaccepting powers in determining the chemical 
reactivity of the ions in aqueous solution. It is expected 
that our co ± s will be helpful to achieve a better under- 
standing of chemical properties of ions in solution and 
further be used in many aspects of solution chemistry 
such as the design of solution-phase reactions according 
to these indices. 
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